import control    #pip install control
import matplotlib.pyplot as plt
import math
import numpy as np

TimeArray = [i*0.001 for i in range(5000)]
K = 2#比例系数
T = 0.5 #时间系数
sys = control.tf([K*T,K],[0,0,1]) #系统
TimeArray, yout = control.step_response(sys,TimeArray)

plt.title("Simulation")  # 括号当中输入标题的名称
plt.plot(TimeArray,yout, label="K=%.f T=%.f" % (K,T))
plt.legend(loc='best')  # 图列位置，可选best，center等
plt.xlabel('Time(s)')
plt.ylabel('Output')
plt.show()


unsuccessful
# 纯微分不能物理实现，也就是仿真的模型的分子的阶次不能比分母的高
# https://www.ilovematlab.cn/thread-33763-1-1.html